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We discuss (2+l)D topological phases on non-orientable spatial surfaces, such as Mobius strip, real 
projective plane and Klein bottle, etc., which are obtained by twisting the parent topological phases 
by their underlying parity symmetries through introducing parity defects. We construct the ground 
states on arbitrary non-orientable closed manifolds and calculate the ground state degeneracy. Such 
degeneracy is shown to be robust against continuous deformation of the underlying manifold. We 
also study the action of the mapping class group on the multiplet of ground states on the Klein 
bottle. The physical properties of the topological states on non-orientable surfaces are deeply related 
to the parity symmetric anyons which do not have a notion of orientation in their statistics. For 
example, the number of ground states on the real projective plane equals the root of the number 
of distinguishable parity symmetric anyons, while the ground state degeneracy on the Klein bottle 
equals the total number of parity symmetric anyons; In deforming the Klein bottle, the Dehn twist 
encodes the topological spins whereas the Y-homeomorphism tells the particle-hole relation of the 
parity symmetric anyons. 
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One of the most striking discoveries in condensed mat¬ 
ter physics in the Ic^Gfew decades has been topolog¬ 
ical states of matterElQ, in which the characterization 
of order goes beyond the paradigm of Landau’s sym¬ 
metry breaking theory. Topological states are gapped 
quantum many-body systems in which the topologi¬ 
cal order is ck^acterized by their fractional excitations 
called anyons&B. These emergent anyons have non-trivial 
braiding statistics, i.e., dragging an anyon once around 
another causes a unitary transformation on the quan¬ 
tum state of the system. In particular, the braiding of 
anyons gives rise to a phase factor in abelian topological 
states. With a non-trivial anyon content A, topologically 
ordered states exhibit long-range entanglemenl^and a 
robust ground state degeneracy (GSD) on torusQs. 

If the microscopic Hamiltonian underlying a topolog¬ 
ically ordered state possesses a symmetry dictated by a 
symmetry group G, such symmetry can interplay with 
the topological order described by A, giving rise 
ferent symmetry enriched topological (SET) state^HiJ. 
If a symmetry A4 in G is on-site, it permutes the anyon 
labels in such a way that the braiding statistics is pre¬ 
served. However, if a symmetry A4 in G involves parity, 
the anyon relabeling changes the braiding statistics by 
a conjugation. In addition, such parity symmetry also 
changes the anyon position. We take the parity to be 
acting in cc-direction in this paper; so the ^-coordinate of 
any anyon is flipped under the parity symmetry. 

One way in which a symmetry can intertwine with the 
topological properties of a state is to twist tker-theory 
by the symmetry through introducing defectalJcS. As¬ 
sociated with each symmetry M in G, we can introduce 
a point-like twist defect from which a branch cut em¬ 
anates. Any anyon passing across the branch cut gets 
transformed or twisted by the symmetry action A4. Un¬ 
like the anyons which are finite-energy excitations, the 
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defects can only be extrinsically imposed by suitably 
modifying the microscopic Hamiltonian. Moreover, the 
defects can braid with projective non-abelian statistics 
and trap exotic zero modes. The study of physical prop¬ 
erties of the defects helps us to gain understanding about 
the topological state itself in the presence of the symme¬ 
try group G. For example, the different possible ways 
by which defects can fuse and associate with each other 
classify the SET states under Glljlij. 

Twisting by parity opens a window for us to study par¬ 
ity symmetric states defined on non-orientable surfaces. 
Similar to on-site symmetries, parity can also be used to 
define defects and to twist topological phases. Though 
the construction of defects for the parity is similar to that 
for the on-site symmetries, the parity defects differ from 
the on-site defects in that they affect the anyon position. 
In addition to the anyon relabeling, the parity defects also 
flip the ^-coordinate for any anyon passing through the 
branch cut. In other words, such branch cut acts like flip¬ 
ping over the surface spatially. Consequently, by suitably 
adding the parity defects, we can turn the underlying sur¬ 
face non-orientable. While topological actioB^igenerally 
fail to be defined on non-orientable surfacefjEil, intro¬ 
ducing parity branch cuts allow us to discuss topological 
phases on non-orientable surfaces. 

A key question which we aim to answer in this pa¬ 
per is what we can learn by putting a parity symmet¬ 
ric topological state on non-orientable surfaces; We may 
gain understanding to the state-it^f by examing how it 
responds to different topologiesQoES. This is partly mo¬ 
tivated by previous studies on the edge states supported 
by (2-|-l)D parity-protected topological phasesEjcZl. In 
these studies, the edge theory, which is a conformaljhfiM 
theory (CFT), is put on a non-orientable spacetzmeEjEEl 
by parity twist, defining an orientifold CFT. Putting the 
edge theory on a non-orientable spacetime allows us to 
infer the trij^lity or non-triviality of the corresponding 
bulk statecSO: the orientifold edge theory suffers from 
quantum anomalies when the bulk is non-trivial whereas 
it is anomaly-free when the bulk is trivial. Here we aim 
to figure out what we can learn about the bulk state by 
putting it on non-orientable spatial surfaces.. 


A. Overview 

In this paper, we focus on (2 -|- 1)D parity symmet¬ 
ric abelian topological states, though the discussion can 
be extended to non-abelian states. Our goal here is to 
study the physical properties which show up when a state 
is put on non-orientable surfaces. In our setup, the non- 
orientable spatial surfaces are generated by introducing 
parity defects. We first formulate these parity defects 
in terms of microscopic degrees of freedom so as to have 
a microscopic understanding of the non-orientable sur¬ 
faces. Then we consider putting a state on a ribbon 
geometry and subsequently twist the theory by parity 
symmetry through introducing a parity branch cut. In 


the parity-twisted theory, we show that the orientation 
of any loop gets flipped when the loop is carried once 
around the ribbon, meaning that such twisted ribbon is 
actually a Mobius strip M. Once the state is put on a 
Mobius strip M, we find that the anyons which are invari¬ 
ant under the parity symmetry show up as the relevant 
anyons under the twist. Interestingly, these parity sym¬ 
metric anyons exhibit a loss of orientation, in the sense 
that their counter-clockwise statistics is identified with 
their clockwise statistics. 

To proceed, we consider putting the state on non- 
orientable closed surfaces. By compactifying the Mobius 
strip M, which emerges as a result of twisting by the 
parity symmetry, we get the real projective plane P. 
While the torus T generates the orientable closed sur¬ 
faces, similarly, the real projective plane P generates the 
non-orientable closed surfaces. For example, gluing two 
real projective planes together gives the Klein bottle K 
whereas gluing three real projective planes together gives 
the Dyck’s surface. Hence we can construct any arbitrary 
non-orientable closed manifold with the real projective 
plane P. It is well known that topological states ha™ 
topology-dependent GSD on orientable closed surfacesOB. 
Here, we construct the ground states and calculate the 
GSD on arbitrary non-orientable closed manifolds. The 
main result for the GSD is summarized in Eqs. and 
(^). The GSD depends on both the topology of the sur¬ 
face and the parity symmetry we used in constructing 
the surface. In particular, the GSD on the real projec¬ 
tive plane P equals the root of the number of distinguish¬ 
able parity symmetric anyons (For the d efinit ion of dis¬ 
tinguishable symmetric anyons, see Sec. |H C| ), whereas 
the GSD on the Klein bottle IK equals the total number 
of parity symmetric anyons. Besides, we also show that 
the calculated GSD is consistent with the Dyck’s theo¬ 
rem, and hence is robust against continuous deformation 
of the non-orientable surface. 

Implications of our main result on the GSD on non- 
orientable closed manifolds can be well described by com¬ 
paring the following two topological states of matter: the 
double-semion state and the Z 2 toric code, where the par¬ 
ity symmetry in each the two states is unique. These two 
states share the same GSD on the torus, GSD(T) = 4, 
and the same total quantum dimension, and hence can¬ 
not be distinguished by the topological entanglement en¬ 
tropy. When put on non-orientable surfaces, however, 
their GSDs now distinguish these two states. For ex¬ 
ample, GSD(P) = 1 for the double-semion state, while 
GSD(P) = 2 for the toric code. Such difference owes 
to the fact that the number of distinguishable parity 
symmetric anyons in the two states are different. Be¬ 
sides, GSD(IK) = 2 for the double-semion state, whereas 
GSD(IK) = 4 for the toric code. Their GSDs are different 
on the Klein bottle IK because they have different num¬ 
ber of parity symmetric anyons. Actually, their GSDs 
are different on any non-orientable closed manifold. 

Finally, we study the physical action of the symme¬ 
try group of the manifold called the mapping class group 
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(MCG). On the torus T, the MCG is generated by the 
modular S-transformation and the Dehn twist. By study¬ 
ing the action of the symmetries on the ground states, one 
can obtain a matrix representation of the MGG. In that 
case, the S-transformation tells the braiding statistics of 
the anyons whereas the Dehn twist tells the exchange 
statistics of the anyons. Motivated by the calculations 
on the torus T, we do the same thing on non-orientable 
closed manifolds. Knowing that the real projective plane 
P has trivial MGG, we consider the first non-trivial case 
which is the Klein bottle K. The MGG of the Klein bottle 
K is generated by a Dehn twist and a Y-homeomorphism. 
Our result for the matrix representation of the MGG is 
shown in Eq. (|^) and ( | 6 ^ . From the matrix representa¬ 
tion of the MGG, we find that the Dehn twist tells the ex¬ 
change statistics of the parity symmetric anyons and the 
Y-homeomorphism acts like particle-hole conjugation. 


B. Organization of the paper 


The content of the paper is organized as follows. In 
Sec. we review the K-matrix theory which serves as a 
mathematical framework in describing (2 -|- I)D abelian 
topological states and we incorporate symmetries into the 
framework. Als o, w e introduce the notion of symmetric 
anyons. In Sec. Ill, we define the defect branch cut as¬ 
sociated to each symmetry operation A4. In addition, 
we show that any anyon getting across the branch cut 
is twisted by the symmetry A4 through examining the 
boundary condition imposed by the branch cut. In Sec. 
0 - We construct the Mobius strip M by suitably adding 
a parity branch cut on a ribbon, showing that the parity 
symmetric anyons become the only relevant anyons as a 
result of twisting by parity symmetry. Then we compact- 
ify it to get the real projective plane P and hence all the 
non-orientable closed surfaces. In each of the cases, we 
construct the ground states and calculate the GSD. Then 
we show that the GSD is robust gainst continuous de¬ 
formation of the manifold. In Sec. M, we study the action 
of the MGG on the ground states of the Klein bottle K 
and obtain a matrix representation for the MGG. Sec. VI 
provides concrete examples for the theoretical discussions 
made. We also illustrate that, in certain scenario, a pair 
of parity defects can be viewed as a deformed genonllSEJ. 
Appendix provides a detailed discussion for the parity 
symmetric abelian states. We show a structural property 
for the set of parity symmetry anyons which turns out to 
be very useful in the derivations. 


II. ABELIAN TOPOLOGICAL STATES 

We start by introducing the field theoretical frame¬ 
work, the K-matrix theory, for abelian topological states 
in (2 -|- 1)D. Then we proceed to incorporate symmetries 
in the K-matrix theory. At the end of this section, we 
introduce the concept of symmetric anyons. 


A. K-matrix Theory 

An abelian topological state is a state which host 
anyons obeying abelian statistics. Any abelian state is 
described by an abelian Ghern-Simons theory known as 
the K-matrix theory, and ^ specified by an invertible 
symmetric integer matrix ATJ. Gonsider an abelian topo¬ 
logical state characterized by AT with Dim K = N, which 
is described by the following (2-|-l)D Lagrangian: 

^Buik = (1) 

where is an TV-component compact C/(l) gauge field. 
It describes a phase of matter which is gapped in the 
bulk. In the presence of a boundary, such abelian Ghern- 
Simons theory induces a bosonic Lagrangian at the edge, 

-CEdge = -^da;4>^Kdt4>, ( 2 ) 

47r 

where 4> is an N-component bosonic field compactified 
by 27rZ^. In Eq. we have omitted the velocity term 
which comes up from the microscopic physics at the edge. 
Such lagrangian describes N branches of gapless edge 
modes. Each positive eigenvalue of AT corresponds to a 
left-moving branch and each negative eigenvalue corre¬ 
sponds to a right-moving branch. 

In the topological state AT, any quasi-particle excita¬ 
tion is labeled by an integer vector I € , in which 

the quasi-particle fusion corresponds to vector addition. 
Local particles are the constituent particles in the micro¬ 
scopic Hamiltonian and they occupy the sublattice ATZ^. 
Identifying any two excitations differed by a local parti¬ 
cle, we define the anyon lattice as 

A = Z^IKZ^. (3) 

The anyon lattice A forms an abelian group with order 
IdetAj under the fusion rule. In the bulk, the Wilson 
operator = expil'^ carries anyon I along a 

path 7 . If 7 possesses open ends in the bulk, the Wilson 
operator generates a gapped particle-hole excitation. On 
the other hand, there is no energy cost if 7 forms a loop, 
or ends at the parity boundary. At the spatial boundary, 
the held operator ' 0 * (x) = e®* annihilates the anyon I 
at the position x of the edge. By braiding and exchanging 
the anyons, we get the mutual- and self-statistics, 

9w = mod 27r, (4) 

01 = K~^l mod 27r/7r. (5) 

Local particles are mutually bosonic in the sense that 
Ow = 0 mod 27r. In addition, any local-particle is either 
self-bosonic with 9i = 0 mod 27r or self-fermionic with 
9i = 0 mod TT. So in Eq. (||), the self-statistics of anyons is 
dehned up to 27r for bosonic systems and it is determined 
only up to TT if there exists fermionic local particles. Let 
V = |detA:|i/2 is the total quantum dimension of the 
state, the S and T matrices are dehned by 

Sw=A^^^'lV and Tw=6wA^^, ( 6 ) 
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where the T matrix is defined up to a sign if there exists 
fermionic local particles. Any abelian state K is modular 
in the sense that S is unitary. If a state is modular, the 
S and T matrices together with the charge conjugation 
matrix Cw = Sjii, where I = —I is the particle gonjuga- 
tion of the anyon I, satisfy the modular relationsH. In the 
K-matrix formalism, all the statistical properties of the 
anyons are encoded in the matrix K. 


B. Formulation of Symmetries 

Here we discuss the formulatio n o f . ay i Uim etry within 
the framework of K-matrix theorytilEBEBEj. Any trans¬ 
formation M is defined by specifying its action on the 
local particle operators. Such action can be rewritten 
as a transformation on the bulk and edge field variables 
in the K-matrix theory. Consider a state K defined on 
upper x-y plane with an edge along the x-axis. Let 
v>^ = (t;0, stii, v'^) for any vector . In general, 

Ma^{x'')M~^ (7) 

M4>{x)M~^ = U'^4>{sx) ttAT^^x, (8) 

where U € GLn(Z) (unimodular), s = ±1 and y S 
Hence, each M is characterized by the data (C/, s, y). 
For example, the data (±7, 1, 0) characterize the trivial 
transform I and the particle-hole conjugation C respec¬ 
tively. While any JVl with s = I acts in an on-site man¬ 
ner, any A4 with s = — 1 involves a reflection from x^ to 
x^ and is called a parity transformation. The operation 
A4 implies a transform on the particle operators, 

MW^M-^ = ( 9 ) 

M'il}\x)M-^ = tp^^sx) ( 10 ) 

Such transformation naturally induces a map on the 
quasi-particles, that is, AIzw : —>■ >->■ Ul. Since 

[/ is a unimodular matrix, its inverse U~^ exists and is 
also unimodular. Hence such mapping is an automor¬ 
phism on the quasi-particle lattice Z^. 

With the understanding of the action of the operations, 
we are ready to define the symmetry of a state. Any 
abelian state K is symmetric under M if the bulk and 
edge actions are invariant under the transformation (|^). 
Such condition translates to the equation 

sUKU'^ = K, (II) 

which is independent on the vector y. If the state K is 
symmetric under M, then M is called a symmetry of 
K. Note that any state K possesses on-site symmetries 
I and C. However, the existence of parity symmetry is 
not guaranteed. If par to symmetry M exists, the signa¬ 
ture for K must vanishEj, meaning that there is an equal 
number of left movers and right movers at the edge. Gen¬ 
erally, parity symmetry, if exists, may not be unique and 
any two parity symmetries differ by an on-site symmetry. 


Besicto, if parity symmetry A4 exists, it must be of even 
ordeiH, or otherwise the K matrix vanishes identically. 
In the paper, unless otherwise specified, we consider par¬ 
ity symmetries of order two. 

If A4 is a symmetry of a state K, then the automor¬ 
phism A4zn on the lattice Z^ descends to a well-defined 
automorphism A4^ on the anyon lattice A, 

Ma ■■ a ; I ^ Ul. (12) 


Since is an automorphism on the anyon lattice A, it 
can be thought of as a relabeling of the anyons. Under 
such relabeling A4a, the S and T matrix transform as 


Suiuv 


Tuiuv 


Sw if 

s = 1 

Stv if 

s = —1 

Tw if 

s = 1 

if 

s = —1 


(13) 


It is known that S and T matrices measure the counter¬ 
clockwise statistics whereas S* and T* matrices measure 
the clockwise statistics. Under a parity symmetry A4, 
orientation of braiding and exchange is flipped due to the 
parity operation. Hence a full parity symmetry transfor¬ 
mation of S and T matrices is a relabeling followed by 
a complex conjugation. Therefore, S and T matrices are 
preserved under any symmetry M.. In other words, any 
symmetry Ai of the topological state K is a symmetry 
of the anyon statistics. 


C. Symmetric Anyons 

In this subsection, we introduce the notion of symmet¬ 
ric anyons, which show up later when we twist the phase 
by the symmetry M. If Ad is a symmetry of the state 
K, we have a group automorphism A4a on the anyon lat¬ 
tice. In the anyon lattice, we have a subgroup consisting 
of symmetric anyons which are invariant under A4a, 

<S = { a € A I A4^(a) = a }. (14) 

While viewed in the integer lattice Z^, any symmetric 
anyon a is only changed by a local particle under the 
relabeling, that is, Ua = a mod KA. Let Ui G Z^ and 
u = {ui, U 2 ,..., um} be a set of generators for S such 
that any a € S can be written as a = aiUi mod KA 
where € Z. Then we can rewrite the subgroup as 


S = uZ^/KZ^, (15) 

where {u)ii = {ui)i is a. N x JL7 matrix with column 
vectors given by the generators^ By a slight abuse of 
notation, we refer, by a generating set, to a matrix of 
column vectors. Albeit all the symmetric anyons can be 
distinguished statistically in A, some of them statistically 
look the same in S. Take any a, a' G S, we define 

a a' if 0aa" = Sa'a" Vo" G 5, (16) 
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choosing the gapping term that couples the two edges, 
we can create a defect branch cut which glues the cut. 
In particular, the gapping term construction for parity 
branch cut gives us a microscopic view on the construc¬ 
tion topological states defined on non-orientable surfaces. 


A. The Defect Branch Cut 


FIG. 1. (a) The state K is cut along the x-axis. The edge 
1 and edge 2 as a whole is considered as a combined edge. 
Physically, the two constituent edges do not talk to each other, 
(b) Without adding any gapping potential coupling the two 
edges along the combined edge, a Wilson line in the bulk can 
end at any position on the two edges. 

where the equal sign for the statistical angles is defined 
up to mod 27r. It gives a equivalent relation which relates 
symmetric anyons having the same statistics in the sub¬ 
group S. Identifying symmetric anyons having the same 
braiding statistics in subgroup S, we have the group of 
distinguishable symmetric anyons in S to be 

5 = 5/^. (17) 

To make it clear, we define the group homomorphism 
fu\s '■ S a uK~^a. In terms of the group 

homomorphism, the equivalent relation can be written 
as a ~ a' if a — a' € Ker/„|s and hence S = iS/Ker/„|s 
is a quotient group. By the fundamental theorem of ho- 
momorphisms, the group S is isomorphic to /u|s[<S] = 
. The number of distinguishable sym¬ 
metric anyons in S is simply given by the number of 
lattice points in such latticeEj. If S saturates S, then all 
the symmetric anyons are statistically distinct, mining 
that the subgroup S itself forms a modular theoryCl. 

III. TWISTING BY A SYMMETRY 

Here we consider cutting and gluing a state K along 
the x-axis. Along the cut, we have gapless modes con¬ 
tributed from the two edges. It is known that by adding 
a gapping term, we can impose a boundary condition on 
the quasi-particle operators at the interface. By suitably 


Imagine cutting a state K along the i-axis (Fig. 0(a)). 
The cut can be physically considered as a combined edge 
along the x-axis. We are going to show that, associated 
to each symmetry M of the state K, we have a gapping 
term ACm for the combined edge that couples the two 
constituent edges by introducing hopping of local parti¬ 
cles with a twist by the symmetry Ai . Denote the lower 
and the upper region as 1 and 2 respectively. For the 
edges 1 and 2 , we have variables (pi and (p 2 respectively, 
which are described by the lagrangians 

Cl = -^d:^(p^{x)Kdtpi{x) (18) 

47r 

C 2 = S^dx4>2{s^)c^^t(p2{sx). (19) 

47r 

Without coupling the fields on the two edges, Wilson 
lines on one side do not talk to those on the other side 
and hence can end at any point on the edges (Fig. 0(b)). 
Let $(a:) = {(pi{x), <p 2 {sx))'^ , the combined edge along 
the x-axis is described by the total lagrangian 

C = Ci+C2 = ^a,$^(x)Ot$(x), (20) 

47r 

where K. = —K 0 sK is the enlarged K-matrix with di¬ 
mension 2N. Let 6 / be an N components vector with the 
only non-zero component at the /-th entry and equals 
one. Take A/ = (e/, since sUKU^ = K, 

AjiCAj = -ejKej + seJU-^KU'^-^ej = 0, ( 21 ) 

for any I, J = 1, 2,..., N. Hence {A/} satisfies the 
Haldane gapping criterion and is hence a possible choice 
of null vectors for /Ce 3. So corresponding to the symmetry 
A4, we can construct a gapping potential 


ACm = —21 cos Af /Cd>(x) — Trefx =^— 2tcos ejK{(pi{x) — U^(p2{sx) — ttK ^x) 

I I 

= (x) + V’f (x)^(MV’f (x)M-i) 


( 22 ) 


which is actually a hopping term for the local particles across the constituent edges of the cut. While hopping from 
the lower edge to the upper edge, ACm twists the local particles by the symmetry M. We will see later that ACm 
also twists aiumns by the symmetry Al, and hence the gapping term ACm defines the defect branch cut for the 
symmetry Adu. 
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Next, we 
orthogonal. 


are going to show that the gapping term ACm is symmetric under the symmetry operation Ad if 17 is 
By using the cosine expression of the gapping term in Eq. , we have the transformation 


MACmM-^ 


—2t cos 

I 

—2t cos 

I 

y^ —2t cos 

I 


ej K M4>2isx)M~^ — ttK~^x) 

e’l KU'^{(j)i{sx) — U'^(j) 2 {x) — TrK~^x) 

seJU~^K{(j)i{x) - U'^(l)2{sx) - ttK~^x) i 


(23) 


where we changed the variable from x to sx in the last 
step. Since U is orthogonal, that the image of the basis 
{ei, 62 ,..., Bn} under U is just the basis itself up to 
some sign changes of the basis vectors. Because sign is 
not important in the argument of cosine, we have 

MACmM-^=ACm- (24) 

So ACm is symmetric under the symmetry Ad, adding 
the defect branch cut does not break the symmetry Ad 
of the whole system. In particular, the matrix U can 
always be chosen as an orthogonal matrix for any parity 
symmetry of order two, hence such parity branch cut is 
always symmetry preserving. 

B. Boundary Condition Across the Branch Cut 

Here, we derive the boundary condition for the Wilson 
operators across the branch cut imposed by the gapping 
term ACm- By coupling the two edges, ACm connects 
the Wilson lines on both sides of the interface. We are 
going to show that the gapping term ACm twists the 
anyons by the symmetry Ad across the branch cut. Since 
the gapping term ACm pins the argument of the cosine 
term to the local minimum of minus cosine, quantum 
mechanically, the possible eigenvalues of the operator, 

G{x) := (l)i{x) — (f> 2 {sx) —'kK~^Xj (25) 

is given by 2'kK~^p, where p is an amzpn carried by a 
Wilson line along the defect branch cutc3. Pick an eigen¬ 
state with fixed p, we have the following equation 

V’r^(sx)V'l(x) = ( 26 ) 

Since {sx)il)\{x)) ^ 0, it describes a particle conden¬ 
sate of —I at-jr^n edge 1 and Ul at sx on edge 2 along 
the interfacald’O. Physically, it should be expected that 
an anyon I getting to edge 1 from below at x can fuse 
with a condensed particle pair to becomes an anyon Ul 
getting away from edge 2 at sx. 

Consider Wilson operators which do not involve any 
Wilson line along branch cut. In such case, p is fixed and 
any physical operator must not change the eigenvalue of 
G{x). More precisely, any physical operator must com¬ 
mute with G{x). It can be easily checked that, without 


a partner on the other side, any Wilson line on one side 
ending on the edge does not commute with G{x). Con¬ 
sider two Wilson lines, one on each side, ending on the 
edges 1 and 2 respectively as indicated in Fig.l. Pulling 
back to the edges, the Wilson operator becomes the field 
operator (xi). For such Wilson operator to 

be physical, we have the condition that 

['(p2\^2)tpi{xi),G{x)] = 0 

/ 2 Sgn(x — SX 2 ) — Ulisgn^x — xi) = 0 

I 2 = Uli and X 2 = sxi. (27) 

Pictorially, while crossing the interface from below, any 
anyon I is twisted by U, at the same time, the anyon 
trajectory goes from x to saH. For any on-site sym¬ 
metry A4, the anyon trajectory is unaffected across the 
branch cut, so the gapping term ACm glues the upper 
side and the lower side directly with a twist U of anyon 
label (Fig. ^a)). So loop orientation is preserved across 
any on-site defect branch cut. The gapping term ACx in¬ 
troduces a trivial gluing and corresponds to a direct hop¬ 
ping term across the interface. The gapping term ACc 
condenses particle pairs which is equivalent to adding a 
SC wire on the interface. For any parity symmetry A4, 
the anyon trajectory is twisted from x to —x across the 
branch cut, hence the gapping term ACm glues the up¬ 
per side and the lower side with a spatial twist in addition 
to the anyonic twist U (Fig. ||(b)). Generally, the change 
in loop orientaion together with the anyon relabelling en¬ 
sure that the anyon statistics is locally well-defined across 
any branch cut. 

The discussion above concerns the boundary condition 
for the Wilson operators crossing the branch cut. Here 
we remark on the boundary condition for implied by 
the gapping term ACm ■ From Eq. (p^ ) , we have 

4)i{x) = M4>2{x)M~^ -I- 2t:K~^p, (28) 

which relates the edge variables (j)i and ^2 at the inter¬ 
face. By differentiating such equation w.r.t. x, we get the 
boundary condition for Ux across the branch cut, that is 

aixix) = Ma2x{x)M~^. (29) 

Hence the branch cut defined by ACm twists by Ad 
from edge 2 to edge 1. In other words, the gapping term 
ACm imposes a twisted boundary condition for the bulk 
gauge field at the defect branch cut. 
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FIG. 2. The gapping term /S.Cm glues the two edges by im¬ 
posing boundary condition on the Wilson operators crossing 
the branch cut. Essentially, it twists anyons by the symmetry 
M across the branch cut. (a) The gapping term for an on-site 
symmetry glues the edges directly with an anyonic twist U. 
(b) The gapping term for a parity symmetry glues the edges 
with a spatial twist together with an anyonic twist U. 


FIG. 3. (a) The state K is prepared in a square shape with 

length scale L. The upper and lower edges are labeled by 1 
and 2 respectively, (b) The twisted ribbon is constructed by 
passing the upper edge through the back and gluing it to the 
lower edge using the gapping term 1\Cm ■ If is on-site, 

the twisted ribbon is a cylinder C; If corresponds to a 

parity symmetry, it is a Mobius strip M 


IV. STATES ON THE NON-ORIENTABLE 
SURFACES 

In this section, we are going to construct topological 
states on non-orientable surfaces. We construct a twisted 
ribbon by adding a defect branch cut on a ribbon geom¬ 
etry along the axial direction and show that the twisted 
ribbon is a Mobius strip if the defect branch cut corre¬ 
sponds to a parity symmetry. On the Mobius strip, we 
see that only the parity symmetric anyons are left. To 
proceed, we define topological states on the real project 
plane by compactifying the Mobius strip. Using the real 
project plane as the generator, we can define topological 
field theories on any non-orientable closed surface. We 
construct the ground states on arbitrary non-orientable 
closed surface and obtain the ground state degeneracy 
(GSD). Using the Dyck’s theorem, we then discuss the 
robustness of such GSD by showing that it is invariant 
under any smooth deformation of the surface. 


A. The Twisted Ribbon 

In this section, we construct the twisted ribbon which 
is a ribbon with a defect branch cut along the axial di¬ 
rection and it gives a Mobius strip if the branch cut cor¬ 
responds to a parity symmetry. We discuss the physical 
Wilson loop operators generally on the twisted ribbon, 
showing that only the symmetric anyons are relevant un¬ 
der the twist. Then we discuss the physical properties of 
the parity symmetric anyons on the Mobius strip. 


1. Physical Wilson Loops Operators 

Consider a state K with a symmetry A4. The state is 
prepared in a square shape with width L in x-direction 
and length L in y-direction (Fig. ||(a)). Edge 1 at y = L 
is passed through the back of the state and glued to the 
edge 2 at y = 0 by using the gapping term A£m ■ We call 
the resulting ribbon a ribbon twisted by the symmetry 


A4. Consider dragging a loop along the non-trivial cycle 
of the twisted ribbon. If ACm is on-site, the orientation 
of the loop is unchanged after a round trip, hence the 

S ted ribbon is topologically the same as a cylinder 
On the other hand, if ACm is a gapping term for 
parity, then the loop orientation is flipped under a round 
trip and therefore the twisted ribbon is a Mobius strip M. 
The two surfaces share the same first homology group, 

i7i(C) = i7i(M) = ( a ) = Z, (30) 

where a is a single cycle going in positive y-direction. 
A loop is either one-sided or two-sided, depending on 
whether its regular neighborhood is an annulus or a 
Mobius strip. The cycle a in the cylinder C is two-sided 
whereas the cycle a in the Mobius strip M is one-sided. 

On the twisted ribbon, any physical Wilson loop op¬ 
erator must have its final anyon the same as its initial 
anyon in the anyon lattice. We are going to show the 
that set of all physical Wilson loop operators is precisely 

W = { IU“ I a G 5 }. (31) 

In other word, only the symmetric anyons are relevant. 
Under the operation IU“ * W^' = IU“+“', which is the 
multiplication of Wilson operators up to a trivial phase, 
W forms a group isomorphic to the subgroup S. 

Notice that all operators in W are physical, so it suf¬ 
fices to show that any physical Wilson loop operator can 
be written as an element in W. Any Wilson operator 
going along a trivial cycle is physical and can be rep¬ 
resented by the trivial element in >V. Any Wilson 
operator going through the cycle a hy k times can be 

written as *■■■* 

In order for such Wilson operator to form a loop physi¬ 
cally, the initial and the final anyons must be the same 
in the anyon lattice. Hence we have the requirement that 
UH = I mod KK which implies [I + U + - + G S. 

Therefore, any physical Wilson loop operator going along 
the cycle a by fc times is also an element in W. To con¬ 
clude, the set of all possible physical Wilson loop opera¬ 
tors is precisely given by the abelian group W. 























2. The Parity Symmetric Anyons 


On the Mobius strip M, the only relevant anyons are 
the anyons invariant under the parity symmetry, which 
we call the parity symmetric anyons. Mathematically, 
we can always find a Mobius strip M on any arbitrary 
non-orientable surface. Hence parity symmetric anyons 
show up on any non-orientable surface. Consequently, 
the topological properties of a topological state defined 
on non-orientable surfaces are deeply related to these par¬ 
ity symmetric anyons. Here we discuss the properties of 
these parity symmetric anyons on M and the structure 
of the subgroup S. 

Firstly, the statistics of the parity symmetric anyons 
exhibits a loss of orientation. Under the parity symme¬ 
try, recall that the S and T matrices are invariant un¬ 
der the anyon relabeling followed by a complex conjuga¬ 
tion which takes into account the change of orientation of 
statistics. Notice that the S and T matrices while they 
are restricted to S are invariant under the relabeling, so 

Saa' = s:,, and Tw=T;V, (32) 

where the equation for the T-matrix is defined up to a 
sign for systems which involve local fermions. Physically, 
it means that the parity symmetric anyons have their 
left hand and right hand statistics identified, exhibiting 
a loss of orientation. By the definition in Eq. (|^) , we get 
a constraint on their mutual-statistics 


daa' = 0 or TT mod 27r. (33) 


Hence the parity symmetric anyons are either mutually 
bosonic or semionic. Similarly, using Eq. (^), we get a 
constraint on their self-statistics 


0 or TT mod 27r, for bosonic system; 
0 or ^ mod tt, otherwise. 


(34) 


Eor a bosonic system, any parity symmetric anyon is ei¬ 
ther a self-boson or a self-fermion. If the microscopic 
hamiltonian involves local fermions, we may have self- 
semionic parity symmetric anyons. More pictorially, 
imagine braiding or spinning a parity symmetric anyon 
at difference points along the one-sided cycle a in the 
Mobius strip M. While it is carried back to its origi¬ 
nal position, the anyon label is unchanged whereas the 
orientation of the braiding or spinning is flipped. Since 
locally we have a well-defined statistics, the initial statis¬ 
tics must agree with the final statistics. So the counter¬ 
clockwise statistics merges with the clockwise statistics 
for the parity symmetric anyons. 

Secondly, there is an explicit formula for the number 
of distinguishable parity symmetric anyons in S. Denote 
[ ] sandwiching a matrix as taking mod 2 to the entries. 
Let tt be a generating set for iS@, the order of S is 

1^1 ^ (35) 


Roughly speaking, each independent column vector in 
[2u^iF~^rt] spans two states. To show this, notice that S 


is isomorphic^ to the lattice [2w^K~^u\'Zi^, it suffices to 
count the number of lattice point in such lattice. Erom 
Eq. (^), the entries of the matrix are in Z 2 . 

We first write such matrix in its Smith normal form R 
which is diagonal and with entries also in Z 2 . Since each 
diagonal entry of R with value 1 mod 2 gives two lattice 
points, the total number of lattice points is Since 

Smith normal form is rank-preserving, Eq. (|^ ) follows. 

Lastly, taking into account the fact that the parity 
symmetry is of order two, as proven in Appendix ^ there 
is an additional structural property for the subgroup S of 
symmetric anyons. If a state K (dim K = N) possesses 
a Z 2 parity symmetry M, then there exists a set v of 
n = N/2 mutually bosonic integer vectors and a set w of 
n' > 0 mutually bosonic integer vectors such that 

u = vLlw (36) 

generates S. The set v is precisely the integer eigenbasis 
of the transformation matrix U with eigenvalue -1-1. We 
call It = r; U w as the generating set for S induced by the 
Z 2 parity symmetry. In term of such generating set, 

K~^v v'^K~^w \ 

K~^v ViFK~^w j 

where K~^v \s a, n x n integer matrix and vo^K~^w is 
an' xn' integer matrix whereas K~^w = 
is anxn' matrix with entries living in ^Z. An immediate 
observation from Eq. ( p^ is that, since sandwiching the 
matrix K~^u by any integer vector must give an inte¬ 
ger, the parity symmetric anyons cannot be self-semionic, 
meaning that they must be either self-bosonic or self- 
fermionic. Another observation from Eq. (H) is that the 
diagonal blocks of the matrix [2u"^K~^u] vanish and the 
entries of the off-diagonal blocks are in Z 2 . Notice that 
the rank of [2u^K~^u] is twice the rank of [2w’^K~^v]. 
From Eq. (^^, the number |iS| of distinguishable parity 
symmetric anyons, like the total number of anyons in a 
Z 2 parity symmetric state, must a square of some inte¬ 
ger. Let V and W be the subgroup generated by v and 
w respectively. Actually, the number of distinguishable 
parity symmetric anyons equals the square of the number 
of anyons in V distinguishable in W^. As a final remark 
on Eq. (^), the column vectors in v (or w) are mutually 
bosonic, so the anyons in the group V (or W) must also 
be mutually bosonic. This statistical property is very 
important later when we construct the ground states on 
non-orientable surfaces. 


u^K-^u = 


B. The Real projective plane 

Here, we are going to construct a topological state de¬ 
fined on the real projective plane P. Generally, to define 
topological states on non-orientable surfaces, we glue the 
surfaces suitably by using the parity branch cuts. Since 
geometrically gluing with a double reflection is a direct 
gluing, to inherit this property for the parity branch cut, 
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FIG. 4. (a) The real projective plane P is constructed by 

shrinking the edge of the Mobius strip M to a point, (b) By 
continuous deformation, the real projective plane P can be 
turned into a cross cap which is a punctured with anti- 
nodal points of the hole identified. The cross cap can be 
viewed as a pair of parity defects on S^. 


we require that twisting twice by parity symmetry A4 
is trivial, namely A4^ = I. Note that the Mobius strip 
has an edge with length 2L. Compactifying the Mobius 
strip M by shrinking the edge to a point (one-point com- 
pactification), we get the real projective plane P, where 
there is no boundary (see Fig. ^(a)). Alternatively, one 
can construct the real projective plane P by starting with 
the square-shaped topological state shown in Fig. ||. The 
method is to glue the top edge to the bottom edge in 
antiparallel sense by using the gapping term IS.Cm for 
parity and to glue the right edge to the left edge in an¬ 
tiparallel sense similarly. The two construction of real 
projective plane are topologically equivalent. By con¬ 
tinuously deforming the manifold, we can turn the real 
projective plane P to a cross cap which is a punctured 
with anti-nodal points of the hole get identified (see Fig. 
^(b)). Such cross cap can be viewed as a pair of parity 
defects connected by a parity branch cut on §^. Let a 
be the one-sided loop going through the cross cap once. 
The first homology group of P is 

iLi(P) = ( a I 2a = 0 ) =Z 2 , (38) 


where a is the generator and 2 a = 0 is the group relation. 
A loop on P is two-sided if it gets through the cross cap 
even number of times. Since 2a = 0, there is no non¬ 
trivial two-sided loop on P. The intersection number of 
the loop a with itself is /(a, a) = 1 mod 2 which is Z 2 
valued since the binary operation / is alternating. Due 
to the intersection of the loop a with itself, we can have a 
non-trivial Wilson operator algebra on the real projective 
plane P. 

To construct the ground states on the real projective 
plane P, we need to figure out the group of physical Wil¬ 
son operators W(P) and then construct a generating set 
for W(P) in which the measuring and raising operators 
can be identified. Since any physical Wilson loop opera¬ 
tor going along a must carry a parity symmetric anyon, 
the set of physical Wilson loop operators is given by 


n;(P) = { I a e 5 }. (39) 

forms an abelian group under the fusion operationEl. For 
any abelian state K wit h parit y symmetry M. of order 
two, recall from section IV A 2 that the set u = v U w 
generates the subgroup S of parity symmetric anyons. 



FIG. 5. (a) The physical Wilson loop operators Wa in the 
group >V(P). (b) By the structural property of the subgroup 
<S, we can construct generators Wa'" and IF“” which can be 
treated as the raising and measuring operators respectively. 


In other words, for any a € S, we can write a = Uy + 
for some ay G V and € W. So the group >V(P) is 
generated by the Wilson operators W^" and FF““, 

W(P) = ( I a, e F, e IT ). (40) 

We have shown that anyons in V are mutually bosonic 
and anyons in W are also mutually bosonic. So we can 
simultaneously diagonalize the operators IT“” and treat 
Wa" as the raising operators (see Fig. ^). We can as well 
construct the ground states another way round but the 
results for the GSD are the same. 

Take the measuring operators IT“™ and the raising op¬ 
erators Wa" , we are ready to construct the ground states 
which corresponds to different fluxes measured by 1T““. 
Let Wi be the i-th vector in the generating set w. From 
the Wilson operator algebra, we have the equation 

IT^W““ (41) 

where vo^K~^ay is the flux vector measured w.r.t. w. 
The measurement of flux vectors w.r.t w defines a group 
homomorphism fyj\v ■ V —>■ ;ay i-T w^K~^ay. 

Take a reference state |0) with W^' |0) = |0). A quantum 
state |a„) is constructed by operating the raising operator 
Wa" on the reference state |0). Two states |a„) and |a(,) 
are identified if they give the same flux. So we have 

|a„)=IT““|0) and (42) 

W^^\ay)=e’-^^<^~"^"\ay), (43) 

where a„ S V/Kei-fiy\v — /i«|y[14]. The first equation 

defines the ground states on the real projective plane P 
while the second equation tells the flux measured by W^' 
for each state. The quantity GSD(P) can be obtained by 
simply counting the number of possible flux vectors in the 
image /iu|y [V] = uF /Z^ . By using the result in 
Eq. (^, it can be shown thafS fw\v[y]'^ — S, therefore 

GSD(P) = (44) 

which is the root of the number of distinguishable parity 

symmetric anyons. Different ground states simply cor¬ 
responds to different fluxes measured by IT““ under the 
flux insertion by W^" ■ Hence GSD(P) counts the number 
of anyons in V distinguishable in W. 
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FIG. 6. The Klein bottle K can be constructed by a connected 
sum of two cross caps. a\ and 02 are the one-sided loops as¬ 
sociated with the first and the second cross caps respectively. 
The Klein bottle K can essentially be viewed as two pairs of 
parity defects on 





FIG. 7. (a) The generators and for the group of 
physical Wilson loop operators on the Klein bottle, (b) By 
the structural property of the subgroup S, the modified gen¬ 
erators IF“ = W^l * 11 ^ 2 “’" and can be treated as the 
raising and measuring operators. 


C. The Klein Bottle 

In the previous section, we have constructed a state 
defined on the real projective plane P. By considering the 
connected sum of two cross caps, we get a state defined on 
the Klein bottle K (Fig. ^). It can be viewed as two pairs 
of parity defects on S^. Alternatively, one can construct 
the Klein bottle K by suitably gluing the square-shaped 
topological state shown in Fig. The method is to glue 
the top edge to the bottom edge in antiparallel sense by 
using the gapping term IS.Cm for parity and to glue the 
right edge to the left edge in parallel sense by using the 
trivial twist A£i. It can be shown that the arrangements 
of parity branch cuts from the two constructions are the 
same. For the sake later of generalization, we stick to 
the picture where the Klein bottle K is a connected sum 
of two cross caps. For each of the two cross caps on the 
Klein bottle K, we have an associated one-sided loop. 
Let «! and 02 be the one-sided loops going through the 
first and the second cross caps respectively. The first 
homology group of the Klein bottle K is 

iLi(K) = ( oi, 0.2 I 2(ai -|- a. 2 ) = 0 ) = Z X Z 2 . (45) 

On the Klein bottle K, there are non-trivial two-sided 
loops. The subgroup of two-sided loops in Hi{K) is gen¬ 
erated by the loops 2ai, 2a2 and /3i = ai + a 2 . The 
intersection number between the loops in Hi(K) can 
be obtained by making use of the bilinear intersection 
I{ai, aj) = Sij mod 2 . 

Now we construct the ground states on K. Similar to 
the previous calculation, we figure out the group >V(K) 
and get a generating set where the measuring and rais¬ 
ing operators can be identified. Any physical Wilson 
operator along a loop in iLi(K) can carry a symmetric 
anyon, so the generating set has IF“^ and where 
a S iS. Since any physical Wilson operator along a 
two-sided loop can carry arbitrary anyon, the generat¬ 
ing set also contains where b € A. 

The operators and carry anyon b counter- 

clockwisely around the first and the second cross caps re¬ 
spectively. The operator carries the anyon b counter- 
clockwisely from the second cross cap to the first cross 
cap and it carries Ub from the first cross cap back to 
the second cross cap. The complete generating set for 
W(K) is Note that 


II /4 = VF-“ + IF|^, where i = 1, 2 and 

(I -I- U)b S S. The generators and can 

be expressed in terms of and , hence they can be 
eliminated from the generating set of yV(IIC). Therefore, 

W(K) = ( a e 5, & G A ). (46) 

Such generating set is not good enough since are not 
mutually commuting though are. Recall that for any 
a G <S, we have a = for some ay € V and Ou, G W. 

Denote IF“ = * W-/™. Since IF“^ * IFg;“”' = IF“, 

the generators IF“^ can be replaced by the operator VF“, 

W(K) = ( w:, I a„ G K, G W, 6 G A ). (47) 

Recall that the anyons in V and W are mutually bosonic 
saparately, so IF“ by such construction are mutually 
commuting. Hence we can simultaneously diagonalize 
FFg^ and treat IF“ as the raising operators (Fig. ^ in 
constructing ground states. 

The ground states corresponds to different fluxes mea¬ 
sured by Let e be the standard basis for Z^. The 
Wilson algebra on the Klein bottle K is given by 

. (48) 

The measurement of flux vectors K~^a defines a group 
homomorphism /e |5 : S —>■ Q^/Z^;a i-T eK~^a. Con¬ 
sider a reference state |0) with IF|(|0) = |0). A state |a) 
is constructed by operating the raising operator on 
the state |0). Two states |a) and \a') are the same if they 
get the same flux vector. So we have 

|a) = W“|0) and (49) 

W|'|a) = (50) 

where the anyon a G 5 /Ker/e |5 ~ The quantity 

GSD(IK) corresponds to the number of flux vectors in the 
image /e| 5 [iS] = K~^u'Z^ . Since e is the standard 

basis for Z^, e can be considered as the identity matrix 
of N dimension. Hence we get /e| 5 [*S] — S. Therefore, 

GSD(IK) = |5|, (51) 

which is the total number of parity symmetric anyons. 
Different ground states corresponds to different fluxes 
measured by the loop operator under the flux in¬ 
sertion by IF“. 
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FIG. 8. Generators for the group of physical Wilson loop operators on q¥ where q is odd. The generators in the first and the 
second column respectively correspond to the generators in the first and the second line in Eq. The figures on the first 

column show the raising operators and the figures on the second column show the corresponding measuring loop operators. 
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FIG. 9. Generators for the group of physical Wilson loop operators on gP where q is even. The generators in the first and the 
second column respectively correspond to the generators in the first and the second line in Eq. The figures on the first 

column show the raising operators and the figures on the second column show the corresponding measuring loop operators. 


D. General Non-orientable Closed Surfaces 

The construction of states on real projective plane 
allow us to discuss states on general non-orientable 
closed manifolds. Mathematically, the set of closed non- 
orientable manifolds up to homeomorphism forms a com¬ 
mutative semigroup under the connected sum The 
semigroup is generated by the real projective plane P, 

{Non-orientable Closed Surfaces} = ( P ), (52) 

where connected sum with P is the same as adding a cross 
cap. The surface P^P is the Klein bottle K whereas the 
surface P^P^P is the Dyck’s surface. With a state on 
real projective plane P, we can generate a state on any 
non-orientable closed surface. Let gP be the connected 
sum of g > 1 cross caps and let ai be the one-sided loop 


associated with the i-th. cross cap. Its hrst homology is 

Hi{q¥) = a, I 2aT = 0 ) = x Z 2 , (53) 

where ax = ai + ■ ■ ■ + aq. Let Pi = ai+ a^+i and FFg^ 
be the operator carrying b counter-clockwisely from the 
i -h 1-th cross cap to the f-th cross cap and Ub from the 
i-th cross cap back to the i -I- 1-th cross cap. Then 

W(gP) =( ..., Wl_^ \aGS, be A). (54) 

Note that 2aT = 0 in Hi{qF), and it requires to be 
trivial. Since the entries of [2u'^K~^u] are in Z 2 , 
commute with the generators . Note that triv¬ 
ially commute with FFg. ’s since there is no intersection 
between ar and PPs. Hence commute with all Wil¬ 
son operators and the group relation is guaranteed. 
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1. Odd Number of Cross Caps 


The construction of ground states for surfaces with odd number of cross caps is a generalization of the construction 
for the surface with one cross cap. Given a non-orientable closed surface with odd number of cross caps, we have 

= 1 = ■ • ■ = 1 = Hence we can replace the generator in Eq. by the operator 

Further rewriting the generating set so as to identify the measuring and the raising operators, we get 


>V(gP) 






IE-, 

IT- 


Uy G V, dyj G 14^, 

bG A 


(55) 


Notice that the generators in the first line commute and the generators in the second line also commute (see Fig. |^). 
So we can simultaneously diagonalize the later generators and treat the former generators as the raising operators. 
Remarkably, non-trivial Wilson algebra appears only in each column and hence the Wilson algebras are decoupled. 
Each of the first (g — l)/2 algebras give |detRr| states while the last algebra gives GSD(P) states. So for q is odd, 


GSD(#) = |detA"|(«-i)/2|.S|i/2^ 


(56) 


The ground states correspond to different fluxes distinguished by the measuring operators shown in second column of 
Fig. g. Take g = 1, the result above reduces to that of the real projective plane P. 


2. Even Number of Cross Caps 


The GSD on manifolds with even number of cross caps is a generalization of the GSD on the surface with two 
cross caps, which is the Klein bottle. Gonsider a surface with even number of cross caps. In this case, we have 

* • • • * instead. So we can replace the operator in Eq. (|5^ ) by the operator 

To figure out the measuring and the raising operators, we further rewrite the generating set as 


W(gP) 


wi 




w 


wk-.^ 

iH- i-Pq-3 


w° 


dy € T, Qjyj 

be A 



(57) 


where IT“ = . Note that the generators in the first line commute and the generators in the second line 

also commute (see Fig. ||). So we can simultaneous diagonalize the generators in the first line and treat the generators 
in the second line as the raising operators. Note that the Wilson algebras are decoupled into columns, where each of 
the first (q — 2)/2 algebras give IdetKl states and the last algebra gives GSD(K) states. So for q is even, we have 

GSD(gP) = |detK|(«-2)/2|5|. (58) 


Each ground state is labeled by the fluxes identified by the measuring operators shown in the second column of Fig. 
If g = 2, the GSD reduces to that of the Klein bottle K. 


From the results in Eqs. (^^ and (^^, we can extract 
the quantum dimension of a cross cap by taking the large 
g limit. Since adding two cross caps generally gives rise 
to an extra IdetiPj amount of independent ground states, 
the quantum dimension of a cross cap is simply given by 
|detK|^/^, which must be an integer and is exactly the 
same as total quantum dimension of the state K. 


E. Robustness of the Ground State Degeneracy 

Here we are going to show the robustness of the GSD 
by checking the consistency of our result with the ge¬ 
ometric Dyck’s theorem. The set of closed surfaces up 
to homeomorphism forms a commutative monoid under 


connected sum ff. The identity element is the sphere, 
and the monoid is generated by torus T and real projec¬ 
tive plane P with a single relation, that is 

{Glosed Surfaces} = ( T, P | P#P#P = P#T ). (59) 

The relation P^K = P^T is known as the Dyck’s theo¬ 
rem. Geometrically, ffT add a handle to the surface with 
the two ends attached on the same side while ffK. add a 
handle to the surface with the two ends attached on the 
opposite sides. In the presence a real projective plane 
P, the surface is non-orientable, so there is no difference 
between ffT and ffK. Let E{g,q) be a closed surface 
given by a connected sum of g tori and g cross caps. The 
Dyck’s theorem is equivalent to the statement that 

S(g,g + 2) = E( 5 -b l,g), (60) 
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for any <? > 1. So the non-orientable surfaces Yi{g,q + 
2) and S(g + l,q) are homeomorphic. We are going to 
prove that their GSD are the same so as to check the 
consistency of our results with the Dyck’s theorem and 
show the robustness of the GSD at the same time. 

To show the consistency, we make use of the following 
observation. The generating set of iJi(S(g,g)) composes 
of the generators of the g tori and the generators of the q 
cross caps. Let pj and (Jj be the standard generators for 
the j-th handle and ai be the generator for the Tth cross 
caps with /3i = ai + a^+i. The generating set for the 
group of physical Wilson operators on E(g, q) are gener¬ 
ated by the union of { Wp ^, } and 

{ ..., } where a € S, A°Notice 

that former generators for the g tori do not intersect with 
the later generators for the q cross caps, so the Wilson 
algebra for the g tori and the q cross caps decouples. So 

GSD(E( 5 , q)) = GSD(I]( 5 , 0))GSD(E(0, q)), (61) 

where it is well known that the factor GSD(E((;,0)) = 
IdetAjs. By making use of the GSD we obtained for 
general non-orientable closed surfaces, we can evaluate 
the other factor GSD(E(0, g)). Consequently, we have 
the following consistency equation 

GSD(E(g, q + 2)) =GSD(E( 5 ,0))GSD(E(0, g))|detiL| 
=GSD(E (5 + l,g)), (62) 

for any g > 1. It shows that the GSD we obtained previ¬ 
ously is consistent with the Dyck’s theorem, i.e., GSD is 
a homomorphism from the monoid {Closed Surfaces} to 
the monoid of positive integers. Physically, the GSD is 
robust against deformation of two cross caps to a genus 
on a non-orientable surface. 


V. THE MAPPING CLASS GROUP OF THE 
KLEIN BOTTLE 

In this section, we are going to discuss the physical 
meaning of the set of large diffeomorphisms, i.e., the map¬ 
ping class group (MCG), of the non-orientable manifolds. 
Notice that the real projective plane P has trivial MCG, 
we focus on the first non-trivial case which is the MCG 
of the Klein bottle. First, we introduce the MCG on the 
Klein bottle and its group action on the group of physi¬ 
cal Wilson loop operators. Then we proceed to its group 
action on the multiplet of ground states. Subsequently, 
we obtain a matrix representation and hence a physical 
interpretation for the MCG. 

A. The Mapping Class Group 

Here we introduce the MCG of the Klein bottle and 
talk about how its group elements act on the physical 
Wilson loop operators in yV(K). Let Aut(X) be the 


group of automorphisms of a closed surface X. Denote 
Auto (A) as the subgroup of automorphisms that are iso¬ 
topic to the identity. The MCG of X is defined as 

MCG(A) = Aut(A)/Auto(A), (63) 

which is the group of isotopy classes of automorphisms of 
X. As proven by Dehn and rediscovered by Lickerish, the 
MCG of any-t^ntable closed surface X is generated by 
Dehn twist^O. In addition, Lickerish has shown that 
the MCG of any non-orientable closed surface X with g > 
1 is generated by Dehn|t«aS|ts and one Y-homeomorphism 
(or the cross cap slide)c3cl. In particular, 

MCG(K) = ( T, Y I TY = YT, T^ = Y^ = / ) (64) 

where the generator T is the Dehn twist and the gener¬ 
ator Y is the Y-homeomorphism. Under the Dehn twist 
T, the Klein bottle is twisted along the two-sided loop 
Pi- On the other hand, under the Y-homeomorphism Y, 
the first cross cap passes through the second cross cap 
and get back to its original position. Essentially, the first 
cross cap has gone through a cycle of a Mobius strip in 
a round trip. 

Note that the Dehn twist T and the Y-homeomorphism 
Y are automorphism for the geometric Klein bottle. 
However, it is not necessarily true for a Klein bottle con¬ 
structed by parity defects. Notice that the configuration 
of parity branch cuts on the Klein bottle is always in¬ 
variant under T. It can be shown that the arrangement 
of parity branch cuts on the Klein bottle is invariant un¬ 
der Y iff. A4^ = I. 

Since the symmetry group MGG(K) acts on the sur¬ 
face K, it induces a group action on the group W(K) of 
physical Wilson loop operators. As shown in Eq. (|4q), 
we can take IU“^ and as the generators for yV(K). 
From Fig. |^, we see that the actions on the generators 
under T and Y transformations, 

TIU“^T-i = W-p, . (65) 

YW:^Y-i = W-;, YW^^Y-^ = Wp^\ (66) 

Since the association of an automorphism in MCG(K) 
with a group action on W(K) is a homomorphism. The 
group relation in MCG(K) are trivially satisfied by the 
induced group action. 

B. The Matrix Representation 

In the previous section, we have seen that elements in 
MCG(K) induce a group action on the group yV(IK) of 
physical Wilson loop operators. Consequently, it defines 
a group action operating on the ground states. In this 
section, we are going to calculate the matrix element of 
the generators T and Y for >V(K). 

Recall that the ground states on the Klein bottle are 
defined in Eq. ( ^^ where we have measuri^ operators 
Wp^ and raising operators IU“. From Eq. (|^), we have 
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FIG. 10. The first two rows show respectively the action of the Dahn twist T on and . An right arrow between two 
figures indicates the action under T. An equal sign between two hgures means a deformation of the loop. The second two rows 
show respectively the action on and under the Y-homeomorphism Y. To be precise, the action under Y is split into 
three steps and the two cross caps are labeled by 1 and 2. In the first two steps, the cross cap 1 is passed through the cross 
cap 2 from left to right. The third step is to drag the cross cap 1 back to its original position. 


[T, = 0, SO we can simultaneously diagonalize the 

operators T and Also, by using Eqs. (^), it can 

be shown that which tells how 

the Dehn twist operates on the raising operators for the 
ground states. Therefore, we can immediately deduces 
the matrix element for the Dehn twist T, 


(a'|T|a) = 


(67) 


where T|0) = e*®'^|0). Hence the matrix representation 
of Dehn twist T is diagonal with the diagonal entries tell 
the topological spins of the parity symmetric anyons. 

Similarly, we can obtain a matrix representation for 
the Y-homeomorphism Y. Note from Eq. that Y 
flips the sign of the anyon label on any Wilson operator. 
Therefore, the Y-homeomorphism Y can be thought of as 
the charge conjugation operator operating on the ground 
states. Hence we get the matrix elements 

{a'|Y|a) = (68) 

where Y|0) = e®®''|0). Hence the matrix representation 
of the Y-homeomorphism Y tells the particle conjugation 
for the parity symmetric anyons. 


VI. EXAMPLES 

In this section, we are going to consider two examples 
of abelian states with Z 2 parity symmetry mentioned in 
Appendix ^ The first one is the fermionic state Kf with 
D = 0 and the second one is the bosonic state Kp. First, 
we obtain explicitly the generating set u = v w for 
the subgroup S of parity symmetric anyons. Then we 
count the number of parity symmetric anyons and the 


number of distinguishable parity symmetric anyons us¬ 
ing Eqs. ( p^ and ( ^ respectively. Next, we calculate 
the ground state degeneracies on arbitrary non-orientable 
closed surfaces. Finally, we write down the matrix repre¬ 
sentation of the generators for the mapping class group 
on the Klein bottle K. 


A. The Fermionic State 

In this subsection, we consider the fermionic state Kf 
with D = 0 and we denote such fermionic state as Kf. 
Generally, the state can possess more than one Z 2 parity 
symmetry, and we consider a particular one here. Let E 
he anxn invertible symmetric integer value matrix. The 
K-matrix K = Kf and the matrix U for the Z 2 parity 
symmetry are respectively given by 

^7=(o-e) ^=(/ 0 )’ 

where the eigenvectors of U with eigenvalue -1-1 are given 
by Vi = (ci Ci)^ where i = 1, 2,..., n. Take any 

a = (oi 02 )^ € S, we have Ua = a -\- KA for some 

integer vector A = (Ai A 2 )^. Such relation can be trans¬ 
lated as a 2 = ai + EAi with A 2 = Ai. Hence, up to 
a local particle, any anyon a G S can be written as 
a = aiiVi mod KA, for some oi £ Z”. From the equa¬ 
tion above, we can see that any parity symmetric anyon 
can be written as a linear combination of the eigenvectors 
Vi's. Hence S is generated by the set u = vUw with 

u = {ui,..., Vn} and w = 0, (70) 

where 0 denotes an empty set. With the generating 

set u, we get the structure of the subgroup S of par¬ 
ity symmetric anyons. By Eq. (^^, one can show that 
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S ~ E~'^1P' j'lP'. The total number of parity symmetric 
anyons is given by the number of lattice points, that is 
|iS| = |deti?|. Observe that the generators in u are mu¬ 
tually bosonic, the matrix K~^u] vanishes. By Eq. 
(^), we get 1 5 1 = 1 . meaning that all the parity sypi,- 
metric anyons are statistically indistinguishable in 5^. 
Since the symmetric anyons in the state Kf are mutu¬ 
ally bosonic, all the physical Wilson operators on the real 
projective plane P shown in Fig. ^ are mutually commut¬ 
ing. Since the Wilson algebra is trivial, there is a unique 
ground state on P. Alternatively, there is no measuring 
operator W®” to distinguish different fluxes, hence 

GSD(P) = 1 (71) 

which is as expected from Eq. On the Klein bottle 

K, the measuring operators are whereas the rais¬ 
ing operators are as shown in Fig. The number 
of different possible fluxes measured by is given by 
the number of flux insertion operators which equals 
\V\ = |5|. Therefore, we have the result 

GSD(K) = |detE|. (72) 

which agrees with Eq. (HI)- Generally, adding two cross 
caps on an unorientable manifold gives extra |detA'| = 
jdetEp ground states. Hence the ground state degener¬ 
acy on a manifold with q cross caps is given by 

GSD(gP) = |detE|«-\ (73) 

The result above suggests that the state K / defined on an 
unorientable manifold is equivalent to the state E defined 
on its orientable double coveiH. Beside, we can consider 
the MCG of K. Since all the parity symmetric anyons 
are self-bosons, we have the result 

(a'|T|a) = , (a'|Y|a) = (74) 

which gives the matrix representation of the generators 
of the mapping class group on the Klein bottle K. 

To understand the parity defects we constructed for 
the state AT/ = E 0 —E, w^i consider a more intuitive 
construction from the genorJlj in a bilayer quantum hall 
state Kbi = E (B E with an on-site Z 2 layer permutation 
symmetry r. The transformation matrix for the on-site 
symmetry r in the state Kbi is exactly the same as the 
matrix U for the parity symmetry of the state Kf . Con¬ 
sider a pair of defects for r connected by an on-site branch 
cut placed along the a;-axis and centered at origin in the 
state Kbi- Such pair of defects for thedayer permuta¬ 
tion symmetry can be viewed as a genonEj. Considering 
flipping over the lower layer about the axis a; = 0, the bi¬ 
layer quantum hall state Kbi becomes the fermionic state 
Kf. It can be checked that the on-site gapping term for 
T turns into the gapping term for parity symmetry by 
sending x to —x for the lower layer edge variables in Eq. 
(p^). Therefore, the resulting pair of parity defects can 
be viewed geometrically as a deformed genon. Note that 


the anyon label of any anyon in a layer going though 
the genon and the deformed genon is always unchanged. 
However any oriented loop going through the genon pre¬ 
serves its orientation whereas any oriented loop going 
through the deformed genon flips its orientation. It can 
also be checked that the GSD contributed by a deformed 
genon in Kf is the same as that contributed by a genon 
in Kbi- In fact, both of them are of quantum dimension 
|detE|. 

To end this example for the fermionic state AT/, we 
consider the particular case where E is an integer with 
E = 2. In such case, the state Kf corresponds to the 
double semion state. Remarkably, in such state, the Z 2 
parity symmetry is unique and is shown in Eq. (|^). In 
the double semion state, we have the group of anyons 
A = {1, si, S 2 , where 1 is the vacuum; si and S 2 
are mutually bosonic semions; ^ is a boson formed by 
the fusion of si and S 2 . The parity symmetry relabels 
the anyons as (1, si, S 2 , 0 ■ 52 j si, ^). Hence the 

subgroup of parity symmetric anyons is 5 = { 1 , which 
is generated by u = r; = {^}. Since all the parity sym¬ 
metric anyons in S have trivial mutual-statistics, they 
cannot be distinguish from one another in S and hence 
S = {[!]} in which ^ get identified with the vacuum 1. 
Gonsequently, we have GSD(P) = I and GSD(K) = 2. 
Generally, on the connected sum of q cross caps, we have 
GSD(gP) = 2«-h 

B. The Bosonic State 

In this subsection, we consider the bosonic state Kb- 
Let A be a n X n non-singular integer value matrix. Such 
topological state may possess more than one Z 2 parity 
symmetry and we pick a particular one. The K-matrix 
K = Kb and the transformation matrix U for the Z 2 
parity symmetry are respectively given by 

Tfb = ( J’t 0 ) and = ( V / ) ’ 

where the eigenvectors of U with eigenvalue +1 are given 
by Vi = (0 Ci)^ where i = 1, 2,..., n. For any a = 
(ai 02 )^ G S, we have Ua = a + KA for some integer 
vector A = (Ai A 2 )^. Such relation can be written as 
Cl = iA(—A 2 ) with Ai = 0. Hence Va G <S, we have 
a = a 2 iVi + {^AA 2 0)^ mod KA. Each integer vector 
from ^AA 2 can give rise to a parity symmetric anyon. In 
other words, each A 2 G Z" with [AA 2 ] = 0 can lead to 
a parity symmetric anyon. Let ([ ] o A)A 2 = [AA 2 ] and 
A = {Ai, A 2 ,..., A„/} be a generating set of Ker([ ] oA). 
We have the generating set u = ?; U w with 

V = {vi,..., Vn} and w = {wi,..., Wn'}, (76) 

where Wi = (^AAi 0)^. With the generating set u, we 
can obtain the structure for the subgroup S. By Eq. 
([Tsl), one can show that S ~ A“^Z"/Z” x [AJZ^ . The 
number of lattice points in such lattice gives the total 
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number of parity symmetric anyons, hence we have |iS| = 
|detA|2^^['''l. In the generating set it = ?; U re, it can 
be checked that the integer vectors in v are mutually 
bosonic and the integer vectors in w are also mutually 
bosonic. However, integer vectors in v can braids non- 
trivially with integer vectors in w. More precisely, the 
matrix [2w;^A"^^ii] = [A^]. By using Eq. (35), we have 

|iS| = Since [ ] o H = [H] o [ ], it can be shown 

that [A] generates Ker[H]. Hence Rk[A] = dim Ker[H] = 
n —Rk[H] where the last equal sign follows from the rank- 
nullity theorem. On the real projective plane P, we treat 
as the raising operators and IE“*’ as the raising 
operators as shown in Fig. The non-trivial Wilson 
algebra leads to a GSD determined by Eq. (|44|), that is 


GSD(P) = 


(77) 


On the Klein bottle K, the measuring operator are 
whereas the raising operators are as shown in Fig. |j. 
By using Eq. , we obtain the GSD 

GSD(K) = |detH|2”-^'^I"^l. (78) 

Since adding two cross caps to an unorientable manifold 
contributes extra |detA'| = |detHp independent ground 
states. For surface constructed by a connect sum of q 
cross caps, we have the general result 

GSD(gP) = |detH|‘?-i2”-’^’'[^l. (79) 

For example, consider the case n = 1 where A reduces 
to an non-zero integer and the bosonic state corresponds 
to a I^A topological state. Finally, we consider a generic 
deformation of the Klein bottle K. Notice that = 
gZTTOj A2 ^ jj0nc0 .^^0 g0t the result 

(a'|T|a) = , (a'|Y|a) = (80) 


which gives the matrix representation of the Dehn twist 
T and the Y-homeomorphism Y. 

To end the discussion for the bosonic state, we consider 
the case where A is an integer with A = 2. The resulting 
state is a Z2 toric code where there is a unique Z2 parity 
symmetry and is given in Eq. ([^. The group of anyons 
in the toric code is M = {1, e, m, ip} where 1 is the 
vacuum; e and m are mutually semionic bosons; '0 is a 
fermionic composite particle of e and m. Under the par¬ 
ity symmetry, the anyons transform as (1, e, m, ip) 1—> 
(1, e, m. Ip), meaning that all of the anyons are parity 
symmetric. The group 5 = {1, e, m, ^jis generated by 
u = vKJw, where v = {m} and w = {e}. Since <S = M, we 
have <S = {[1], [e], [m], [ip]} where there is no identifi¬ 
cation between the symmetric anyons. Consequently, we 
have GSD(P) = 2 and GSD(K) = 4. Generally, on the 
connected sum of q cross caps, GSD((7P) = 2'^. There¬ 
fore, on any given non-orientable closed manifold, the 
GSD of the Z2 toric code is always double that of the 
double-semion state. 


VII. SUMMARY AND DISCUSSION 


In this paper, we explore the physical properties of 
topological states put on non-orientable surfaces. We in¬ 
troduce, for any given parity symmetry, the gapping term 
for the parity branch cut, which glues two edges of a state 
with a spatial flip. We construct a ribbon by gluing the 
opposite edges of a square-shaped topological state using 
the gapping term for parity. After the surgery, we find 
that carrying an oriented loop along the cycle of the rib¬ 
bon flips the orientation of the loop, meaning that the 
ribbon is actually a Mobius strip. On such twisted rib¬ 
bon, only those anyons which are symmetric under the 
given parity symmetry are left. These parity symmetric 
anyons, which have their left hand statistics the same as 
their right hand statistics, determines the physical prop¬ 
erties of the state when put on non-orientable surfaces 
constructed by suitable arrangement of the given par¬ 
ity branch cut. In other words, by studying the proper¬ 
ties of the state on non-orientable surfaces, we can learn 
about the physical properties of these parity symmetric 
anyons. For example, the GSD on the real projective 
plane P equals which is the root of the number 

of distinguishable parity symmetric anyons; the GSD on 
the Klein bottle K equals the total number |5| of parity 
symmetric anyons. Generally, adding two cross caps on a 
non-orientable surface leads to an extra IdetAl indepen¬ 
dent ground states, which is exactly the number of inde¬ 
pendent ground states arise in adding a handle. Gonse- 
quently, such GSD is robust against smooth deformation 
of two corss caps to a handle on top of a non-orientable 
surface. Besides, we study the action of the MCG on 
the ground states of the Klein bottle. We hnd that the 
Dehn twist T encodes the topological spins of the parity 
symmetric anyons whereas the Y-homeomorphism tells 
the particle-hole relation of the parity symmetric anyons. 
Finally, we work on the fermion state Kj and the bosonic 
state Kb as examples to illustrate our results. 

While topological actions cannot be defined purely on 
non-orientable surfaces in a coordinate independent way, 
introducing parity branch cuts opens a window for us 
to study topological states on non-orientable surfaces. 
One particularly interesting scenario is the case where 
all anyons in A are parity symmetric, e.g., the toric code 
model. In such case, the parity acts only on the anyon po¬ 
sition without changing the anyon label. In addition, the 
twisted boundary condition for the bulk gaugerfijeld be¬ 
comes continuous up to a physically trivial partEll. Since 
the state K is modular, we have M = 5 ~ 5. Hence 
the GSD on a connected sum of q cross caps equals 
IdetAT^'^^. The consequence of requiring the gauge field 
to be strictly continuous is discussed in a previous 
related worko. In their work, they quantize the BF the¬ 
ory on non-orientable surfaces with the consideration of 
smooth differential forms and densities. Under such re¬ 
quirement, the GSD of toric code reduces from 4 to 2 on 
the Klein bottle K. However, in the context of topologi¬ 
cal phases, even if all anyons in A are parity symmetric. 
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there is no reason in requiring to be strictly continuous 
across the branch cut. 

In discussing the physical Wilson loop operators on the 
surfaces decorated with parity branch cuts, we require 
that any Wilson operator must have its finial anyon the 
same as its initial anyon and any Wilson loop is living in 
the first homology group of the surface. While the former 
condition guarantees that the Wilson loop operators do 
not create any gapped excitations, the later condition 
is imposed for Wilson loops on orientable surfaceH^. 
Generally on any surface, regardless of its orientability, 
the relevant loop group is the first homology group. To 
justify this, it suffices to show that the line integral of 
the gauge field along any homologically trivial loop is 
zero. Note that if a loop is homologically trivial, it is 
the boundary of some 2-chain which is a formal sum of 
oriented 2-simplexes. Since the flux through each of the 
2 -simplex is zero, the flux through the 2-chain is also zero. 
Hence, by Stokes’ theorem, the line integral vanishes. 

There are some open questions along the direction of 
this work. One open question is the derivation of the 
edge theory on non-orientable surfaces with non vanish¬ 
ing boundary, e.g., the Mobius strip M, which possesses 
a single edge. It would be interesting to know how the 


loss of orientation in the bulk affects the spectrum at 
the edge. Another interesting direction is the calcula¬ 
tion of topological entanglement entropy and entangle¬ 
ment spectrum on a non-orientable subsystem, say, the 
Mobius strip M. It is known that the topological en¬ 
tanglement entropy on a disc shaped subsystem is given 
by lii^p, which is independent on the choice of ground 
statcflO. Also there is a one-to-one correspondence be¬ 
tween the physical edge spectrum of the disc shaped sub¬ 
system and low-lying spectrum of its entanglement 
hamiltoniarHcj. An immediate question is whether the 
same thing hold true on a Mobius strip shaped subsys¬ 
tem, which shows up when non-orientable surfaces are 
suitably bipartited. Last but not least, a future direction 
is to generalize the discussion to put topological states 
on higher dimensional non-orientable spatial manifolds. 
By stacking up layers of (2 -|- 1)D topological states and 
suitably introducing coupling betiEem layers, we can ob¬ 
tain (3 -I- 1)D topological statea£3i£3. By using similar 
techniques, we can stack up the parity defects to obtain 
a reflection defect in (3 -I- 1)D. Alternatively, one can 
construct such reflection defect by coupling adjacent lay¬ 
ers with a parity twist. Such reflection defect can be 
used to define layers-constructed topological states on 
non-orientable spatial manifold. 


Appendix A: Abelian States with Z2 parity symmetry 


Here, we review the general data of the K-matrix theory for abelian topological states with a Z 2 parity symmetry. 
Then we talk about a structural property of the subgroup S of parity symmetric anyons. Mathematically, Ad is a Z 2 
parity symmetry of a state K if it is a symmetry of the state K with s = — 1 and C/^ = I. In other words, AT is a 
state with a Z 2 parity symmetry if there exists a unimodular U such that 

-UKU'^ = K and = I. (Al) 


The equations above serves as the constraints for the matrices K and U. If K and f/ is a solution to the above 


equations, then K' and U' is a gauge equivalent solution if 
U' = GUG~^. Up to gauge equivalence, the general solutior 


exists a un imodular G such that K' = GKG^ and 
I to Eq. (M) is given by 


/ 0 

A 

B 

B \ 


/ 

-I 

0 
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0 

0 

I 

0 ) 


(A2) 


where A, B, G, D, E are integer value matrices with D is anti-symmetric and E is symmetric. Note that A is a rif, x nf, 
matrix whereas B and G are both Ub x Uf matrices. Also matrices D and E are both of dimension Uf x nj. Hence 
Dim K = N = 2n where n = Ub + Uf. The state with n/ = 0 is called the bosonic state Kb whereas the state with 
rife = 0 is called the fermionic state Kf. Any local particle must be a boson in Kb whereas the system may include 
local fermions 'm Kf. 

We proceed to talk about the structure of the subgroup S of parity symmetric anyons. More precisely, we are going 
to show that if a state K possesses a Z 2 parity symmetry Ad, then there exists a set v oi n = N/2 mutually bosonic 
integer vectors and a set w of n' > 0 mutually bosonic integer vectors such that u = v U w generates S. The set v is 
precisely the integer eigenbasis of U with eigenvalue -|-1. The set u = u U w is called the generating set for S induced 
by the Z 2 parity symmetry. 

First, we show that u = v U w generates the subgroup S of parity symmetric anyons for some set w of integer 
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vectors. Induced by the symmetry transformation U, we have a natural basis for the free abelian group 


Vb it — 



Vf if 


I 


0 

0 

e-if 

^if 


ih — 



Xf it — 


I 0 

0 

^if 

V 0 ^ 


(A3) 


where ibjj = I, 2,..., n^/j.Let v and x be the collection of first two and last two species of basis vectors respectively, 
where we have n basis vectors in both v and x. By such construction, v is the eigenbasis for U with eigenvalue +1. 
Let Uo be a generating set of the subgroup S. Since u is a subset of 5, the union v U Uo is also a generating set. 
Note that we can expand any vector in Uo in the basis v U x. Since r; is a subset of the generating set v U Uo, we 
can eliminate the v components of the basis vectors in Uo and denote the resultant set as w. Hence we obtain the 
generating set u = v U w for <S, where any vector in w can be expanded in the basis x, that is w C x1^. We denote 
the number of vectors in such w as n' > 0. The subgroup generated by v is denoted as H whereas the 

subgroup generated by w is denoted asW = wZ" jKlJ^. 

Next, we show that symmetric anyons in Y are mutually bosonic and the anyons in W are also mutually bosonic. 
For any a G S, we have Ua = a + KA for some A G Z^. From the defining properties of Z 2 parity symmetry in Eq. 
(Al), we get [/^A = A. Take any a, a' G S, since H^A = A, we have K^KA' = 0. By expanding the expression 


^Ua', we get a formula for the statistics of symmetric anyons under the Z 2 parity symmetry, 


a^K-'^a' 


= -^ia 


^A' 


JT 


A). 


(A4) 


Since U'^A = A, A is an eigenvector of with eigenvalue equals + 1 . Adopting the basis naturally induced by the 
parity symmetry as shown in Eq. ([A2|), for any symmetric anyon a, we can write a and A in column vectors as 



A = 



(AS) 


where the components abi, 062 , Ab G Z”*" and the components o/i, 0 / 2 , A/ € Z”^. In particular, if o € E, we have 
Obi =0, 0/1 = 0/2 and Ab = A/ = 0. Similarly, if o S W, we have O 62 = 0/2 = 0. By making use of the explicit form 
of K, U, a, A in Eqs. (A2) and (AS), the defining equation Ua = a + KA translates to 


Obi — A 2BA/) 


and ^/i — ^/2 — —Xb — (LI -t- A)A/, 


(A 6 ) 


where CA/ = 0. Given the matrix K, the equations above relates the vectors o and A. Substitute the relations above 
into the formula in Eq. (A4), we get the useful formula for the statistics of the parity symmetric anyons. 


K — —-(o^^A'f, + a'^Xb) + [x^EX'f — (o/^A/ + a'f^Xf)^ + -^A^CA/ + X'^CXf + A/(ZI + Il’^)A/^ 

= A A AA/ — ( 0 / 2 A/ + o/^A/)^, (A7) 

where we made use of the facts that CA/ = CA/ =0 and D is anti-symmetric in the second equal sign. In the second 
line, the object inside the big round bracket is an integer which does not contribute to the mutual-statistics. Take 
any 0 , 0 ' G V, then Ab = A(, = A/ = A/ = 0, hence we have a'^K~^a' = 0. Similarly, take any a, a' G W, then 


Ob 2 = o (,2 = 0 , hence we have o^AT "o' 
anyons in W are also mutually bosonic. 


G Z. Therefore, the symmetric anyons in V are mutually bosonic and the 
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